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We present near- and mid- infrared observations on the shock-cloud interac- 
tion region in the northern part of the supernova remnant HB21, performed with 
the InfraRed Camera (IRC) aboard AKARI satellite and the Wide InfraRed 
o ,1 Camera (WIRC) at the Palomar 5 m telescope. The IRC 7 /iin (S7), 11 /im 

(Sll), and 15 /im (L15) band images and the WIRC H 2 v = 1 -> S(l) 2.12 

! um image show similar shock-cloud interaction features. We chose three repre- 

sentative regions, and analyzed their IRC emissions through comparison with H 2 
line emissions of several shock models. The IRC colors are well explained by the 
thermal admixture model of H 2 gas — whose infinitesimal H 2 column density has 

\ a power-law relation with the temperature T, dN ~ T~ b dT — with n(H 2 ) ~ 10 3 

cm" 3 , b ~ 3, and iV(H 2 ; T > 100 K) ~ 3 x 10 20 cm" 2 . The derived b value 
may be understood by a bow shock picture, whose shape is cycloidal (cuspy) 
rather than paraboloidal. However, this picture raises another issue that the 
bow shocks must reside within ~ 0.01 pc size-scale, smaller than the theoreti- 
cally expected. Instead, we conjectured a shocked clumpy interstellar medium 
.{ph \ picture, which may avoid the size-scale issue while explaining the similar model 

parameters. The observed H 2 -u = 1 — >■ S(l) intensities are a factor of ~ 17 — 33 
greater than the prediction from the power-law admixture model. This excess 
may be attributed to either an extra component of hot H 2 gas or to the effects 
of collisions with hydrogen atoms, omitted in our power-law admixture model, 
both of which would increase the population in the v — 1 level of H 2 . 
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Subject headings: ISM: individual (SNR G89.0+4.7) — (ISM:) supernova rem- 
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Introduction 



Shock-cloud interactions in supernova remnants (SNRs) have been studied for their wide 
astrophysical importance, such as the followings. Firstly, we can study the physics and chem- 
istry of molecular shocks (jDraine et al.lll983t iHollenbach &: McKed Il989l ). to be used as a 
basis for interpreting astrophysical phenomena. Secondly, the evolu tion of shocked molecular 
clouds is closely related to the star formation occ urring in them (jElmegreen Sz Elmegreen 



19781 : IVrba et all Il987l : iRevnoso & Manguml l2001f ). Lastly, the evolution of the SNRs in 



Rho k Petre 


1998: 


Shinn et al. 


2007) 



addition, SNR shocks may have a simpler configuration in their interaction region than do 
shocks in outflows, since the latter can be generated by a diverse range of bullets, clumps, 



and/or winds (cf. Fig. 11 in IHollenbach et al.lll989l ). Hence, SNRs provide geometrically 



simpler conditions to study the shock-cloud interaction. 



The shock-cloud interaction m anifests itself through radiative means (jDraine fc McKee 
1993 IWardle fc Yusef-Zadej] 120021): high-velocity wings on low-J CO rotational lines (e.g. 



Koo fc Moon 



1997 



Koo et al 



the far-infrared (e.g. IWatson et al. 



200 11); high- J CO rotational l ines and OH rotational lines in 



1985; Melnick et al. 



987h: radio maser l ines of OH (e.g. 



Frail et al.lll994j ); and infrared emission lines of H2 (e.g. iNeufeld et al.l 120071 ). Among these, 
the infrared emission lines of H 2 are useful in studying the excitation condition o f shocks, 
since H 2 is the most abundant molecule in interstellar clouds (ISnow fc McCallll2006l ) and the 
lines, originating from quadrupole transition, are optically thin and closely spaced through 
~ 1 — 30 /im. In addition, high resolution imaging is easier for H 2 infrared emission lines, 
than for radio emission lines. 

Infrared H 2 emission lines have been frequently observed in SNRs interacting with 
nearby clouds. One noticeable common result is that the H 2 level-population dia gram shows 
an ankle-like curve (cf. F ig . ED) in the range of energy levels, — 25, 000 K (IQliva et al. 



2001 



1990; iRichter et al.lll995al Jbl; iRho et al.ll200ll : INeufeld fc Yuanll2008l : see Fig. 7a in 



Rho et al. 



for the ankle-like curve), which has been usually interpreted as a population of two- 
temperature (~ 10 2 and ~ 10 3 K) H 2 gas in Local Thermodynamic Equilibrium (LTE). 
Such a population ca nnot be explained by a single planar shock model (e.g. Table 2 of 
Wilgenbus et al.ll2000l ). Hence, several mechanisms have been proposed i n order to explain 



this level population distribution, such as a partially dissociative J-shock (IBrand et al.lll988 
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Burton et al. 1989; Moorhouse et a 



stationary shock (IChieze et al 



1998 



Cesarskv et al.lll999h . 



1991 J) , a bow shock (jSmith et al.lll99ll ). and a non- 



However, the interpretation of the H 2 population diagram is still under debate; for in- 
stance, even for t he same target SNR , IC443, all the above three models were pre ferred by 
different authors (IRichter et al.lll995al ; ICesarsky et al.lll999l ; iNeufeld fc Yuanll2008l ). One re- 
striction that hinders a clear understanding of the H2 population diagrams, and subsequently 
of the shock-cloud interaction, is that the observations hitherto performed were not able to 
resolve the interaction features sufficiently to distinguish their structures, such as planar 
shocks, bow shocks, shocked clumps etc. Our target SNR, G89.0+4.7 (conventionally known 
as HB21), provides an excellent example to study the shock-cloud interaction, avoiding this 
restriction, because of its large angular size and proximity. 



HB21 is a large (~ 120' x 90'), middle-aged 



2006 




Bvun et a 


2006) SNR, located 


Tatematsu et al. 


1990|; 


Bvun et al. 


2006) 



1987 



5000 - 7000 yr, lLazendic fe Slane 
of ~ 0.8 - 1.7 kpc Leahy 



Bvun et al.ll2006l ). It has been suspected of interacting with a molec- 



ular cloud, because of its deformed shell-like shape in the radio and the e x istence of nearby 



giant molecular clouds (jErkes &: Dickell Il969l ; iHuang fc Thaddeusl Il986l ; iTatematsu et al. 
1990|). Such an interaction was confirmed with the detection of broad CO emission lines 



near the edge and the center of the remnant (IKoo et al.l l200ll ; iByun et al.l 120061 ) . This de- 
tection also supported the suggestion that the clo ud evaporation might be re sponsible for 
the enhanced thermal X-rays in the central region ( Leahy fc Aschenbach 19961 ). 



We have obtained infrared images of two specific regions of the remnant, where the 
shock-cloud interaction is under way. The near- and mid-infrared band images (~ 2 — 27 
jum), obtained from two infrared cameras respectively onboard the AKARI satellite and the 
Palomar 5 m Hale telescope, show several shock-cloud interaction features. We derive their 
brightness ratios and analyze them on the basis of the H 2 emission model from shocked gas. 
We find that their infrared properties are well described by shocked H 2 gas with a power-law 
distribution of temperatures. 



2. Observations 



We observed two specific regions ("Cloud N" and "Cloud S" in Fig. |2J), where it is 
known that slow shocks (< 20 km s" 1 ) propagate into clouds of density n(H 2 ) ~ 10 3 cm" 



(IKoo et al.ll200ll ). with two different instruments: the InfraRed Camera (IRC. lOnaka et al. 
2007 ) aboard the AKARI satellite and the Wide-field InfraRed Camera (WIRC. Iwilson et al. 
20031 ) at the Palomar 5 m telescope. In this paper, we present the Cloud N data; the Cloud 
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S data will be considered in a future paper. The details of the observations and the data 
reduction on the IRC and WIRC data are described separately below. 



2.1. AKARI IRC observations 



AKARI was designed for both imaging and spectroscopy in the infrar ed ([Murakami et al. 



20071 ). It has two scientific instruments; one, the InfraRed Came ra fIRC.IOnaka et al.ll2007D . 
covers 2-30 /im, and the other, the Far Infrared Surveyor (FIS. lKawada et al.ll2007l ). covers 
50-200 /iin. AKARI has imaged HB21 several times in either pointing or scanning modes 
with the IRC and FIS. In this paper, we analyze the IRC pointing observation data. 

The IRC pointing observations for Cloud N were performed at 2006 December 5-6 
toward (RA, Dec) = (20:47:35.20, +51:12:07.00) in J2000. IRC comprises three channels 
(NIR, MIR-S, and MIR-L), and each of these has three broad-band filters for imaging. 
Among t hese, we employed six filters, two in each channel, for the observations (IRC02 
mode, see lOnaka et al.ll2007l ). The wavelength coverage and the imaging resolution (r) are 
listed in Table [T], together with the pixel sizes of the detector array in each channel. 

All observation al data were initially processed with the IRC Imaging Pipeline (v. 20070104 
Lorente et al.l 120071 ). In this process, several instrumental corrections were made to the 
dithered images, with the corrected images then coadded. The instrumental corrections 
incl ude a dark correctio n, a bad-pixel masking, a distortion correction, and flat-fielding 



sec 



Lorente et al.l 120071 . for detail). We used the "self-dark," a dark image measured for 
each individual observation, rather than the "super-dark," a general dark-image periodically 
measured during the mission. Astrometric information was appended to the coadded im- 
ages by matching the positions of field point-sources with those of 2MASS catalog sources 
( ISkrutskie et al.ll2006l ): the matching tolerance was 1.5 pixels. We used the most recent con- 
version factor, which converts the instrumental flux (ADU) to physical flux (Jy), as noted 
on the AKARI /IRC Data Reduction Support Page0 at 2007 Dec 20. The systematic errors 
(~ 2 — 5%) of the conversion factor are included in the error estimation. 



For compatibility between the images from different bands, we equalized the pixel size 
and the spatial resolution. The images were resampled to have the same pixel-size (1" x 1"), 
employing the public software SWARP^I, and convolved with different Gaussian kernels to 
have the same spatial-resolution (~ 7.43"). Then, to enhance diffuse features, we removed 



^ttp : //www. ir . isas . jaxa. jp/AKARI/Observation/DataReduct ion/ IRC/ConversionFactor_071220 .html 



http : //terapix . iap . f r/ 



point sources with DAOPHOT package (IStetsonl 119871 ) in IRAF; saturated point-sources 
were simply masked out. The final images of Cloud N are displayed in Figure [31 



2.2. Palomar WIRC H 2 observations 

We also carried out near-infrared imaging observations of Cloud N, centered at (RA, 
Dec) = (20:47:36.45, +51:11:42.41) in J2000, using the 2.12 /mi narrow-band filter for the H 2 
v = 1 —>■ S(l) transition with WIRC on the Palomar 5 m Hale telescope on 2005 August 
29. WIRC is equipped with a Rockwell Science Hawaii II HgCdTe 2K infrared focal plane 
array, covering an 8.7' x 8.7' field of view with a ~ 0.25" x 0.25" pixel size. We obtained 
50 dithered images of 30 sec exposure. For the basic data reduction, we subtracted a dark 
and sky background from each individual dithered frame and then flat-fielded. We then 
combined the dithered frames to make the final image. 

The astrometry solution was obtained by matching the positions of 14 field point-sources 
with those of 2MASS catalog sources. The positions were matched within ~ 0.1", which is 
smaller than the 2MASS systematic root-mean-square uncertainty of 0.15". We calibrated 
the image by comparing the magnitude of the 14 field point-sources with the corresponding 
K s magnitude from the 2MASS catalog. The correlation coefficient and the ratio between the 
two magnitudes was 0.9930 and 1.047 ±0.011, respectively. The systematic error (~ 12%) of 
the calibration factor was included in the error estimation. The point sources were removed 
with DAOPHOT as done for the IRC images, and the full-width-at-half-maximum (FWHM) 
of the resultant point-spread-function is ~ 1.2". The final image for Cloud N is displayed in 
Figure [31 



3. Results 



The final images f rom the IRC and WIRC are displayed in Figure [3j together with 1420 



Koo et al. 



Koo et al. 



MHz radio continuum flTatematsu et alJll990h and 12 CO J = 2 -> 1 230.583 GHz ( 
2001) images for reference. The peak positions of the Nl and N2 clouds, where 
(120011 ) observed 12 CO J = 2 — > 1 broad molecular lines, are indicated as crosses on all 
images. 
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3.1. Morphology 

The IRC images show different features from band to band. Bands N3 and N4 are 
dominated by point sources, although N4 also shows a faint diffuse feature near the N2 
cloud. Bands S7, Sll, and L15 all show similar diffuse features over the observed field. 
Several filamentary feature are commonly seen: the large V-shaped feature around the Nl 
cloud; the short-and-horizontal features northward from the Nl cloud; the horizontal feature 
around the N2 cloud; the vertical feature southwestward from the N2 cloud. The last two 
features are also seen in the N4 band. Bands S7, Sll, and L15 also manifest a clumpy feature 
at the west-end of the N2 cloud, which accompanies a cycloidal diffuse feature around it 
(cf. Fig. H]). The L24 band shows rather extended and diffuse features, which are different 
from the relatively sharp features seen in the S7, Sll, and L15 bands; moreover, the bright 
regions in the L24 band do not overlap with the Nl and N2 clouds. 

The WIRC H 2 v — 1 — > S(l) image shows filamentary features, which are very similar 
with those seen in the IRC S7, Sll, L15 bands; for example, the west side of the V-shaped 
feature around the Nl cloud, and the horizontal and vertical features near the N2 cloud. 
This similarity is more recognizable when we compare the WIRC image with the RGB color 
image made with S7 (blue), Sll (green), and L15 (red) (see Fig. [5]). In this RGB image, 
the diffuse features commonly seen in the IRC S7, Sll, L15 bands are more prominent; 
especially, the clumpy feature near the west-end of the N2 cloud is more identifiable. Also, 
the RGB image shows that the infrared "color" varies from red to blue across the image. 
The filamentary features seen in the WIRC image well follow the bluish features seen in the 
RGB image. 

To quantify the similarity of the diffuse features commonly seen in the IRC S7, Sll, L15 
images and the WIRC H2 v = 1 — > S(l) image, we measured the correlation of the features. 
We sampled an area as outlined in Fig. [5] with a yellow polygon excluding yellow circles with 
a red slash, and made scatter plots between the IRC and the WIRC images (Fig. [6]). Each 
point of the scatter plots corresponds to the mean intensity of an 8" x 8" region. As upper- 
panels of Figure El show, the correlation of S7-to-Sll (0.94) and Sll-to-L15 (0.95) are very 
high, while that of S7-L15 (0.84) is relatively low, but still high. The correlations of the IRC 
images to the WIRC image vary from 0.63 to 0.78 (lower-panels of Fig. [6]); the S7 image 
has the highest correlation, while the L15 image has the lowest. This is consistent with the 
above description that the filamentary features seen in the H 2 v = 1 — > S(l) image well 
follow the bluish features seen in the RGB image. The positive y-intercepts indicate that 
the spatially-correlated diffuse features seen in the IRC images are put on some background. 
We removed these backgrounds during the intensity measurement of the diffuse feature, by 
subtracting the intensity of nearby regions (cf. § 13.21 and dashed-line boxes in Fig. 
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Some of the features seen in the IRC and WIRC images contain geometrical relationships 
with the 12 CO J = 2 — > 1 and 1420 MHz radio continuum images. The V-shaped feature 
overlaps with the U-shaped Nl cloud seen in the 12 CO J = 2 — » 1 image, and its apex is 
located near the Nl cloud peak position. Westward from this apex, there is a bright region in 
1420 MHz radio continuum. The horizontal feature around the N2 cloud well matches with 
the N2 cloud seen in the 12 CO J = 2 — > 1 image, and the 1420 MHz radio image also shows 
a diffuse feature there. The clumpy feature at the west end of the N2 cloud also shows an 
isolated clumpy feature in the 12 CO J = 2 — > 1 image, although it is rather vague because 
of the low spatial resolution of the CO image; no counter part of this feature is seen in the 
1420 MHz image. Considering the northward propagation of the remnant (cf. Fig. and 
the geometrical connections with the features seen in the CO and radio continuum images, 
the V-shaped feature around the Nl cloud and the horizontal feature around the N2 cloud 
appear to represent a bow shock and a blast wave, respectively. Similarly, the clumpy feature 
at the west end of the N2 cloud seems to be a shocked clump, considering the propagation 
direction of the remnant, together with the cycloidal diffuse feature surrounding the clumpy 
feature (cf. Fig. EJ). 



3.2. Quantitative Infrared Characteristics of the Shocked Gas 



For the analysis of the shock-cloud interaction features, identified in the previous section, 
we first chose three representative regions for further examination. As the RGB image 
(Fig. |3j) shows, the interaction features have different colors from region to region, which 
means that their infrared characteristics vary. The three regions we chose are designated as 



follows: 



Nlwake," "N2front," 



and "N2clump" (see Fig. ED . 



Then, we measured the intensity of these three regions in the IRC images by subtracting 
the background emission from nearby areas. Two background regions were carefully selected 
for each chosen region, excluding any diffuse filamentary features (see Fig. [5]). To avoid any 
possible contamination from point sources, we also excluded these. The measured intensities 
and the derived color values (S7/S11 and S11/L15) are listed in Table [2j 

The intensity in the WIRC H 2 w = 1 ^0 S(l) image was also determined by subtract- 
ing the backgrou nd emission from nearby areas. Additionally, it was extinction-corrected. 
Lee et al.l (120011 ) estimated the foreground hydrogen nuclei column density, iV(H)=iV(H 
l)+2iV(H2), from the X-ray absorption toward the central region of HB21, and found iV(H) = 
(3.5 ± 0.4) x 10 21 cm" 2 . This JV(H) corresponds to A? 

12// , ™, — 0-22 ma g in the case of the 
"Milky Way, R v = 3.1" curve JWeingartner & Draine!l200ll ; lDraine!l2003h . The tabulated H 2 
v — 1 — ► S(l) intensities in Table [2] have been corrected to compensate for this amount of 
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extinction. 

Table [2] shows that the Nlwake is the brightest in the L15 band, while the N2clump 
and N2front are the brightest in the Sll band. The S7 intensity of the Nlwake is below the 
3-cr limit (cf. Fig. [3]); we indicate its intensity as a 90% upper confidence limit. The color- 
color diagram of S11/L15 vs. S7/S11 (Fig. [7]) reveals that each region has distinctive color 
properties: the N2front is the bluest, while the Nlwake is the reddest. This color property 
is also well-recognizable from the RGB map (Fig. The H 2 d = 1 ^ S(l) intensity is 
the brightest at the N2front. The H 2 v = 1 — > S(l) intensity of the Nlwake is below the 
3-cr limit, again, as in S7 band; its intensity was also indicated as a 90% upper confidence 
limit. From the colors and H 2 i; = 1^0S(l) intensity of the selected regions, we deduce 
that the bluer the color is, the higher the H 2 v = 1 — > S(l) intensity is. 



4. Radiation Source of the Shock-Cloud Interaction Features Observed in the 

AKARI IRC Bands 

Since IRC images are broad-band images, they may contain contributions from both 
continuum and line emission. Thus, to understand the observed infrared color ratios, we 
should first identify the radiation source which generates the shock-cloud interaction features 
observed in IRC S7, Sll, and L15 bands. We consider that the emission lines of H 2 are the 
main source of the interaction features, based on the following arguments. 

Firstly, the features seen in the IRC S7, Sll, and L15 band images are similar to those 
seen in the Palomar WIRC H 2 v = 1 — > S(l) image; the correlation is about 0.63 — 0.78 
(cf. §3.11 and Fig. E]). This strongly suggests that H 2 emission is responsible for at least 
some of the fe atures seen in the mid -infrared IRC bands. Secondly, a recent result supports 



this assertion. iNeufeld et al.l (120071 ) observed the shock-cloud interaction regions in several 
SNRs, and detected rotational lines of H 2 and HD, together with fine-structure lines of 
several atoms and ions. From their emission line maps, it was found that the emission lines 
show distinctive spatial distributions, categorized into five distinct groups. Among these 
groups, only the "lines of S and H 2 (J up > 2)" group spans the emission in the IRC S7, 
Sll, L15 bands (~ 5 — 20 /xm), which suggests that only the emission lines in this group 
can generate the highly- correlated shock-cloud interaction features in the three IRC images 
(cf. upper-panels of Fig. E]). Among the emission lines belonging to this group, only H 2 
lines fall wit hin the three IRC bands th emselves. Finally, the assertion is also theoretically 



supported. Kaufman fc Neufeldl (119961 ) numerically simulated slow C-type planar shocks 
(v s = 10 — 40 km s _1 ), propagating into fully molecular gas of n(H 2 )= 10 4 — 10 6 cm -3 . 



According to their results, for a case similar to Cloud N (v s = 20 km s l , n(H 2 )= 10 4 cm 6 1 
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H2 lines are the dominant emission lines over the 5 — 20 /im spectral range. 



We could consider other possible sources for the emission, such as fine structure ionic 
lines, thermal dust co ntinuum, Polycycli c Aromatic Hydrocarbons (PAHs) bands, and syn- 
chrotron radiation (cf. iReach et al.ll2006l ). However, these do not seem likely, as we discuss 
below. 

Within the wavelength coverage of the IRC S7, Sll, L15 bands, there are two ionic lines, 
[Ne II] 12.8 [im. and [Ne III] 15. 5 /im, whose emitter (Ne + , Ne 2+ ) have bee n observed in the 



shocked regions of SNRs (e.g. iRaymond et al.lll997l ; iNeufeld et al.l 120071 ) . Although these 



lines cannot generate spatially-correlated features in the three IRC band images, they may 
contribute to the shock-cloud interaction features seen in Sll and L15 images. However, 
this does not seem likely because of their poor spatial-correlation with H 2 emissions, well 
correlated with the shock-cloud interaction features seen in the IRC bands (cf. lower panels of 
Fig- E]). [Ne II] and [Ne III] show negative spatial-correlations with H2( Ji m > 2) in three SNR s 
out of four, and the one exception shows a low correlation of ~ 0.4 (INeufeld et al.l 120071 ). 
lower than the correlations of Sll and L15toH 2 v = 1 ^ S(l), 0.7 5 and 0.63 (cf. Fig. El). 
Addi tionally, a shock of > 80 km s -1 is required for Ne to be ionized (IHollenbach fc McKee 
19891 ). whose velocity is far higher than the velocity derived from CO observations, ~ 25 
km s _1 (iKoo et al.ll200ll ). Hence, the contribution of [Ne II] and [Ne III] to the shock-cloud 
interaction features seen in the three IRC bands would be negligible. 



In C-shocks, the grain temperature is below ~50 K (IDraine et al.lll983l ). thus negligible 
thermal dust continuum is generated over the 5 — 20 /im range. Besides, below ~ 50 K, ther- 
mal dust continuum increases toward > 20 /im; hence, if significant thermal dust continuum 
were generated by the shocks, the IRC L24 image would also show the same shock-cloud 
interaction features as in the S7, Sll, and L15 bands, considering the not-that-low sensitivity 
of the L24 ban d. (The L 24 sensit ivity in the IRC02 observation mode is nearly half of the 
L15 sensitivity, lOnaka et al.ll2007l .) This is not the case. Therefore, thermal dust continuum 
is not likely to generate the interaction features seen in the IRC S7, Sll, and L15 images. 



PAHs, thought to be formed by grain shattering in shocks (I Jones et al.lll996l ). show 
strong and broad band feat ures at 3.3, 6.2, 7.7, 8.6, 11.3, and 12.7 /im. However, they have 
not been observed in shocks (Ivan Dishoeckll2004l ; lTielenal2008h . This invisibility is attributed 
to the property that PAHs are heated slowly a nd cool efficie ntly, hence generate little excess 
emission above the background PAH emission (ITielensI 120081 ) . Moreover, C-sho cks of 25 
km s _1 are not strong enough to sputter refractory grain cores, like PAHs (IDraine et al. 



19831 ). Thus, we reject PAHs as a source for the shock-cloud interaction features. 



Synchrotron radiation is unlikely, as well. IKoo et al.l (120011 ) derived a spectral index of 



-10- 



a ~ -0.28 (fluxcx v~ a ) with intensity I (y = 325 MHz) ~ 0.05 K ~ 160 Jy sr" 1 at Cloud N. 
If we extrapolate this intensity up to ~ 10 fim, this corresponds to ~ 7 Jy sr -1 . This is far 
less than the observed intensity of ~ 10 5 Jy sr -1 (see Table |2j). 



5. Comparison to the H 2 Infrared Emission of Shock Models 

Following the arguments in $U we now compare the observed infrared characteristics 
of the shock-cloud interaction features, to that predicted for H2 emission lines from various 
shock models. 



5.1. C-Shock: Isothermal H 2 Gas 

In the Nl a nd N2 clouds, s low C-type shocks (y s < 20 km s~ x ) were observed from the 



CO observation ( Koo et al. Also, for a planar C-type shock, it is known t hat the shock- 



heate d gas can be approximated as an isothermal and isobaric slab of gas (INeufeld et al. 



2006) , in view of the Ho excita tion diagrams predicted for such shocks (e.g. lKaufman fc Neufeld 
19961 ; IWilgenbus et al.ll2000l ). Hence, we first calculate the expected IRC colors from the 
emission lines of isothermal H 2 gas. 



For the model calculation, we follow INeufeld fc Yuanf s (120081 ) assumptions. They an- 
alyzed shock-cloud interaction features seen in Spitzer IRAC images based on H 2 infrared 
emission lines. We assume that the molecular gas consists of H 2 and He only and that the 
density of He, n(He), is proportional to the density of H 2 , n(H 2 ), with n(He)=0.2n(H 2 ). 
Then, we calculate the H 2 level populations in statistical equilibrium, only including the col- 
lisional excitation through H 2 -H 2 and H 2 -He collisions. The ortho-t o-para ratio (OPR) of H ? , 



is assumed to be 3.0. Collisional excitation rates are obtained from lLe Bourlot et al.l (119991 ) 



and other molecular data fo r H 2 were obtained from the database provided by a simulation 
code, CLOUDY (v07.02.01. iFerland et al.lll998l ). The effects of extinction are also included 



when modeling the IRC colors, with iV(H)= 3.5 x 10 21 cm 2 , as in §3.21 

The modeled IRC colors from isothermal H 2 gas in statistical equilibrium are plotted as 
open circles (o) in Figure [7J together with the observed IRC colors of the selected regions. 
Since pure rotational levels of H 2 have their shorter wavelength transitions at higher energy 
levels, the trajectory of the expected IRC colors moves from the lower-left corner to the 
upper-right corner as the temperature increases. Points for the same n(H 2 ) approach to the 
LTE value as n(H 2 ) increases. 



Figure [7] shows that isothermal H 2 gas cannot explain the observed IRC colors with any 



- 11 - 



combination of density and temperature, although the Nlwake may possibly originate from 
isothermal H 2 gas at T~ 300 K, since its S7/S11 value is an upper limit. We also vary the 



OPR fr om 0.5 to 5, since th e OPR is expected to 



clouds (IDalgarno et all Il973l: iFlower fc Watj 



3e different from 3.0 in the interstellar 



shocked gas (ITimmermann 



1998; 



1984 



Wilgenbus et al. 



Lacy et al.l Il994l ) and possibly even in 



2000h . However, the observed IRC colors 



are not reproduced. The predicted IRC colors at the same temperature vary according to 
the adopted OPR; however, the locus of predicted IRC colors are not much different from 
the OPR=3.0 shown in Figure El 

This disagreement is consistent with the well-known fact that the H2 level population 
seen in shock-cloud interaction regions shows an ankle-like curve (cf. fjl]). The critical 
den sity of an H? line tran sition increases as the energy level of the upper state increases 
(cf. iLe Bourlot et al.lll999l ); hence, an isothermal H 2 gas can only produce either a straight 
line (LTE) or a knee-like curve (non-LTE) in the population diagram (cf. Fig. [1]), neither of 
which are observed. 

The disagreement mentioned above is understandable for the N2clump, since it may 
possess a shock structure similar to a bow shock, in which the shock velocity is non-uniform. 
However, understanding this disagreement between the model and the observations is not 
straightforward for the N2front, since morphologically it looks like a planar shock from its 
elongated, uniform appearance, together with the northward propagation of the remnant (see 
Fig. E]). This suggests that another process is needed to understand the IRC colors of the 
N2front under the planar C-shock model. One possible explanation is the Wardle instability, 
which develops in C-shocks w hen an unstab le balance between the magnetic force and the 
ion-neutral drag force breaks (jWardkjll990l ). However, it has been s hown that this process 
results in a negligible cu rvature in the H 2 level population diagram (INeufeld fc Stondll997l ; 
Mac Low fc Smith! 119971 ). contrary to the data. 



5.2. C-Shock: Power-law Distribution of H 2 Gas Temperature 



The observed IRC colors are not reproduced by the planar C-shock model, in which 
the shocked gas is isothermal (cf. §5.11) . Hence, we model an admixture of H 2 gas with a 
temperature distribution. We assume a powe r-law distribution for th e column density based 
on the temperature, dN = aT~ b dT; recently, INeufeld fc Yuan! (120081 ) found that this model 
well describes the infrared characteristics of the shock-cloud interaction in the SNR IC443, 
where H 2 emission lines are dominant. alN is the infinitesimal H 2 column density within 
the temperature range (T,T + dT), a is a constant related to the total H 2 column density, 
iV(H 2 ), and b is the power-law index. We adopt a temperature range for the integration as 
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(T m in, T max ) = (100 K, 4000 K). T min is set to be lower than lNeufeld &: Yuanf s value of 300 K, 
since the AKARI IRC b ands (S7, Sll, L15) are sensitive to lower temperature H 2 gas than 
the Spitzer IRAC bands iNeufeld fc Yuan! (120081 ) used, which only extends to 8 fim. T mnx is 



set a s 4000 K, since H 2 gas rapidly dissociates above this temperature (e.g. 
20021 ). With these variables, a can be written as follows 



Le Bourlot et al. 



iV(H 2 ;T > 100K)(6- 1) 



T 1 



b 

mm 



max 



X 



, where iV(H 2 ; T > 100 K) is a total column density of molecular hydrogen warmer than 100 
K. The effect of extinction is also included, as for the isothermal H 2 g as case. 

The modeled IRC colors for this case are plotted as filled circles (•) in Figure [71 As is 
evident, the observed IRC color ratios can now be reproduced with a suitable combination 
of n(H 2 ) and b. Interestingly, both the N2clump and N2front have a similar b value, ~ 3, 
with n(H 2 ) of ~ 10 2 — 10 3 cm -3 . The column densities are found to be iV(H 2 ;T > 100 K) 
~ 3 x 10 20 cm -2 . The derived parameters from the power-law admixture model are listed in 
Table [3j For each set of parameters, the detail contributions of H 2 line emission to IRC bands 
are listed in Table |U The "Weight" column of Table H] lists the weighting factor for each 
line in the IRC band contribution. For example, the Sll band intensity can be calculated 
as follows. 

WH 2 ) 0.610 J[H 2 5(2)] + 0.921 J[H 2 5(3)] 



MJy sr 1 10 4 erg s 1 cm 2 sr 1 

As Table H] shows, the pure-rotational H 2 emission lines are dominant in all IRC bands. 



(2) 



The derived parameters are similar to those previously determined towards several 
SNRs, where interaction with nearby molecular clouds is occurring. The derived n(H 2 ) at 

l arge V elocity 
fepOlh . From 



^q'^xIO 3 cm 3 , is similar to the value, derived from 



3.1 x 1 3 cm 3 bvlKoo et al. 



the N2front, n(H 2 ) = (l 

Gradient analysis of CO data for HB21, n(H 2 
their Spitzer IRAC observation towards the SNR IC443, INeufeld fc Yuan! (120081 ) found that 
the IRAC color ratios were well explained with a range of p ower-law index b, 3.0-6.0. Our b- 
values (~ 3) fall into near the lower end of INeufeld fc Yuanl 's range. Also, the column density 



we derive, A r (H 2 ; T > 100 K) ~ 3 x 10 20 cm 2 , is similar to iV(H 2 ) towards shock-cloud inter- 
action regions in four other SNRs (W44, W28, 3C391, and IC443), iV(H 2 )= (2.8-8.9) x 10 



20 



cm 2 , which were deter mined from a two-tem perature LTE fitting of pure-rotational H 2 spec- 



tra with varying OPRs (INeufeld et al.l 120071 ). From these three parameters derived (n(H 2 ), 
iV(H 2 ; T > 100 K), and b), we also determined the model predictions for the H 2 v = 1 — > 
S(l) intensities, and list these in TableEJ However, we note that they are a factor of ~ 17—33 
smaller than those observed (see Table [2]). We discuss these results further in $6j 
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5.3. Partially Dissociative J-shocks 



The partially dissociative J-shock model was pr oposed by iBrand et al.l (119881 ) to ex- 



plain the ankle-like curve of H2 level population (see iBurton et al.l Il989l for details on the 
model). In this model, it i s assumed that H 2 gas survives the J-shock jump (< 25 km s^ 1 , 
Hollenbach &: McKedll980l ). and cools behind the shock front. The postshock temperature 
covers T < 3 x 10 4 K, for which t he maximum corresponds to a shock velocity v s h ~ 25 



km s _1 (IHollenbach fc McKedll979l ). The level populations for the surviving H2 gas behind 
the shock are determined by the cooling function, which depends on the temperature and 
the density of the postshock gas. The semi-analytical model only includes H 2 cooling (line 
and dissociational) and CO line cooling, and has two free parameters to adjust: the pressure 
(P) and the fractional CO abundance, X co =iV(CO)/iV(H 2 ). The expected IRC colors at 
different (P, X co ) are shown as open circles (o) in Figure [9) Extinction is also included in 
modeling the IRC colors, as in §5.1[ adopting iV(H) = 3.5 x 10 21 cm -2 . We yary X cn from 



to 10 and varied P from 10 cm K up to 10 cm K. In comparison, [Burton fc Haas 



( 119971 ) obtained the best fit to mid-infrared H 2 observations toward Orion Molecular Cloud 
Peak 1 for P = 8 x 10 10 cm" 3 K. 

Figure [9] shows that none of these model parameters falls within the observed color 
range for any of the three regions in HB21. The closest correspondence is for the N2front, 
with the parameters (P, X co ) = (~ 10 11 cm -3 K, ~ 10~ 6 — 10~ 5 ); this pressure is higher than 
that expected in SNRs, as we now show. The shock velocity needs to be less than 25 km s~ l 
for th e H 2 to survive. It is also known that n(H 2 )~ 10 3 cm -3 for the N2 cloud (IKoo et al. 
200ll ) . Hence, the postshock pressure (~ pv 2 ) of the remnant would be ~ 10 8 cm" 3 K , about 
1000 times smaller than the derived value, ~ 10 11 cm -3 K, above. iMoorhouse et al.l (ll99ll ) 



and IChevalierl (119991 ) showed that the postshock pressure can be higher than the ambient 
pressure, when the radiative shell of the remnant collides with dense molecular clumps; 
however, it is only a factor of ~20 higher. There is also the possibility that the postshock 
H 2 gas does not cool as low as a few hundred K. In this case, the lower temperature H 2 gas 
would be absent, causing the IRC colors to move towards the upper-right direction in the 
color-color diagram (Fig. |7j) to bluer colors. Thus, this could not explain the discrepancy 
between the model and the observations. Overall, a partially dissociative J-shock does not 
seem to be a suitable model to explain the observed IRC colors. 



5.4. Non-stationary C-Shocks 



The non-stationary C-sh ock model describes the temporal evolution of a shock-cloud 
interaction after the collision. IChieze et al.l (119981 ) showed that the non-stationary C-shock 
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model possesses both J- and C-shock charact eristics around an evo lutionary time of ~ 1, 000 
yr, before the steady state has been achieved. ICesarsky et al.l (119991 ) showed that the modeled 
H 2 level populations can be applicable to observations of the shock-cloud interaction region 
in the SNR IC443. We obtain the H 2 line intensity of five pure rotational lines, S(l)-S(5), 
modeled for a shock velocity v s h = 25 km s _1 , a preshock hydrogen nuclei density n H =10 4 
cm -3 , where n H =n(H)+2n(H 2 ), and a preshock magnetic field B = 100/xG at the evolution 



ary ti mes of (1.25, 1.50, 2.00, 4.00) xlO 3 yr, from the results of iFlower fc Pineau des Forets 
( 119991 ) : this presho ck density n H is a little higher than that expected for Cloud N of HB21, 
n(H 2 ) < 10 3 cm -3 (IKoo et al.ll200ll ). The expected IRC colors for these models are overplot- 
ted as filled circles (•) on Figure [9) As the time elapses, the steady state C-shock is reached, 
and the expected color approaches that of isothermal gas; this is consistent with the result 
that the shocked gas be hind a C-shock can be treated as an isothermal and isobaric slab of 
gas JNeufeld et alfeood ). 



Figure [9] shows that none of the expected colors overlaps with the observed IRC colors. 
However, when considering the age of HB21, ~ 7,000 yr, the elapsed time after the shock- 
cloud collision may be less th an ~ 1, 000 yr. Besides, the mod el parameters predict the 



postshock n H to be ~ 10 5 cm 3 ( Flower fc Pineau des Forets! 1999h . which i s 100 times higher 



than derived from the shocked CO gas, n(H 2 ) ~ 10 3 cm 3 (jKoo et al. 2001 ). Hence, it would 



be worth examining the IRC colors at earlier times (< 1, 000 yr) and with a lower preshock 
n(H 2 ) density (< 10 4 c m" 3 ). One way des cribing the non-stationary shock at such early 
times was developed by lLesaffre et al.l (120041 ) . However, predictions for the H 2 line emission 
from this model are not yet available, so this conjecture cannot be tested. 



6. Discussion 

As shown in £j5l only the power-law admixture model can explain the observed IRC colors 
at the N2clump and the N2front; the other three models — the planar C-shock, the partially 
dissociative shock, and the non-stationary C-shock — can not simply reproduce the observed 
colors. Hence, here we discuss the derived parameters for the power-law admixture model, 
based on two pictures of the shock-cloud interaction, in a bow and in a clumpy medium. 
Additionally, we here comment on the adopted OPR value. Since no information on the 
actual OPR is available at the moment, we assumed OPR=3.0; thus, the three parameters 
derived can be changed according to the adopted OPR. 
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6.1. Nature of Molecular Shocks Seen in the Mid- Infrared 



6.1.1. Bow Shock Picture 



Recently, iNeufeld fc Yuan! (120081 ) showed that values they derived for the power-law 
index b, 3.0-6.0, in their power-law admixture model can be explained by paraboloidal bow 
shocks (i.e. a geometrical combination of planar C-shocks). They derived a power-law index 
b ~ 3.8 for a paraboloidal bow shock, when the H2 survives the shock (i.e. T < 4,000 K). 
They argued that if some bow shocks have a maximum shock- velocity which is too small to 
fully cover the temperature range (< 4, 000 K), then the spatially averaged value of b will 
be higher than 3.8. In this model, higher-T H2 gas originates near the apexes of the bow 
shocks, while lower-T H 2 gas originates from the wakes of the bow shocks, since the shock 
velocity (the normal to the bow shock surface) decreases from the apex to the wake. 

Our derived b- values (~ 3), somewhat lower than 3.8, can be understood within the 
bow-shock picture as follows. If the bow shock has a cycloidal shape (cf. Fig. H]), squashed 
enough along the axis of symmetry, the relative size of the slow shock-velocity portion of 
the bow will be decreased in comparison to the paraboloidal case. This will decrease the 
amount of the lower-T H 2 gas, leading to a lower power-law index b. The cycloidal feature 
surrounding the N2clump may imply this possibility, although the feature seems to locate 
behind a compact clump rather than ahead of it as the bow shock should be (see Fig. [5]). The 
N2front does not show any bow-shock-like features in the IRC and WIRC images (cf. Fig. [5]). 
In this case, the derived b value can be understood in the bow shock picture if the bow shocks 
are unresolved. However, this leads to a problem associated with the size-scale for the bow 
shock. 



The distance to H B21 is known to be ~ 0.8 — 1.7 kpc (iLeahyl 119871 ; iTatematsu et al. 



19901 ; iByun et al.ll2006l ) and the imaging resolution of the Palomar WIRC H 2 image is ~ 1.2" 
(cf. §2.2p . This angular size means that the bow shock must be hidden within a length scale 
of ~ 0.01 pc ~ 3 x 10 16 cm. This length is comparable to the width of a single planar 
C-sho ck, propagating in t o a molecular cloud of preshock density n H ~ ^(H 2 ) 



10 3 cm -3 

^e.g. iDraine et al.lll983l ; iTimmermannl Il998l ; IWilgenbus et al.ll2000l ). Since a bow shock is 
a geometrical combination of planar shocks, this scale is too small to hide the bow shock. 
In other words, the observed color ratios and the uniform-and-elongated appearance of the 
N2front suggest that all the H 2 gas, with 100 < T < 4, 000 K, must be mixed within a length 
scale ~ 0.01 pc, however, a bow shock, with a preshock density n H ~ 10 3 cm -3 , is unable to 
mix H 2 gas of such temperatures within such a length scale. 



The bow shock picture has an additional difficulty in supplying t he preshock H 2 ga s. In 
the outflow case, for which the bow shock picture was first proposed ( iSmith et al.lll99ll ). the 
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preshock H2 gas is supplied from the interstellar medium (ISM), which spreads far upstream; 
hence, the preshock H 2 gas may be continuously supplied to the bow shock front. On the 
other hand, in the SNR case, the preshock H 2 gas is mainly supplied from the swept-up 
molecular gas, which only spreads a finite distance upstream. Hence, the amount of preshock 
H2 gas may be insufficient to develop the steady-state bow shock, in the SNR environments. 



6.1.2. Shocked Clumpy ISM Picture 



The derived model parameters for the power-law admixture model may be understood 
under a shocked clumpy ISM picture. When a clump is swept up by a planar shock, a 
deformation of the clump is expected. Indeed, from numerical simulations, it was shown 
that such a defor mation bring abo u t a c ycloid al tail whe n a sm ooth-edge clump is swept 
up. For example, iNakamura et al.l (120061 ) and IShin et al.l (120081 ) showed the formation of 
a cycloidal tail for a non-magnetized and magnetized clump, respectively. Their collision 
conditions are the same; the Mach number of the shock and the density contrast of the cloud 
relative to the intercloud medium are both 10. As is evident, the shapes of these tails are 
very similar to the cycloidal feature seen in the N2clump (see Fig. [5]). 

In view of this similarity, we conjecture the following picture. If the conditions for the 
shock-cloud interactions at the N2front and the N2clump are similar except for the clump 
size, the resultant excitation of the H2 will also be similar. Under this assumption, this 
is quantified by the parameters n(H.2) and b of the power-law admixture model. Then, 
two observational results for the N2clump and the N2front that are not clearly explained 
in the bow shock picture are now simultaneously explainable. First, the similar power-law 
index b ~ 3 results from the similar collisional conditions. Second, the different appearances 
(filamentary vs. clumpy) results from the physical size of the shocked clump. In this picture, 
the size of the clumpy ISM at the N2front must be less than ~ 0.01 pc ~ 3 x 10 16 cm, to 
be unresolved in the WIRC H 2 u = 1 ^ S(l) ima ge (Fig. [5 |). Sug gestively, this size is 
comparable to that of a shocked clump, proposed by I Chevalier! (119991 ) to explain the high 
velocity molecular gas of low column density, observed in the SNR IC443. 

However, this picture also has a issue related with the very linear shape of the N2front 
seen in the H^v = 1 — ► S(l) image. If numerous unresolved-clumps exist in the N2front, the 
overall appearance may be wriggly si nce the shock would propa gate further at the less dense 
regions, as shown in the simulation of iPatnaude &: Fesenl (120051 ). One possible explanation is 
that the size of clump is so small (<< 0.01 pc) that the expected wriggle is also unresolved 
in the H2 v = 1 — » S(l) image. However, it is uncertain whether the cycloidal tail would 
be formed in such a small scale. 



-17- 



6.2. H 2 v = 1 -> S(l) intensity 



The observed mid-infrared IRC colors are well explained by a power-law admixture 
model of H2 gas temperature (cf. Fig. [7]). The dominant emission lines in the IRC S7, Sll, 
L15 bands are five pure-rotational lines, v = — > S(l)-S(5). In order to examine whether 
this model can also explain the shock seen in the near-infrared, we compare the predicted H 2 
v = 1 — > S(l) intensity from the derived parameters — n(H 2 ), b, and iV(H 2 ; T > 100 K) — to 
the observed intensity. The observed and predicted H 2 w = 1 ^ S(l) intensities are listed 
in Table [2] and [3], respectively. As is evident, the observed intensities are a factor of ~ 33 and 
~ 17 larger than the predicted intensities, for the N2clump and the N2front, respectively. 

One simple explanation for this disagreement is the existence of additional H2 gas, whose 
temperature and density are both high, but whose column density is low enough to have 
negligible effect on the mid-infrared line intensities. To compensate for the deficiency of the 
H 2 v = 1 — > S(l) intensity requires N(v — 1, J = 3) ~ 10 1 
if there is additional H 2 gas present, in LTE with T ~ 2000 K and iV(H 2 ) 
both the mid- infrared IRC colors and the near-infrared H 2 d = 1 ^ S(l) intensity can 
be understood. A compact, unresolved shocked cloud is a candidate for producing such H 2 
emission. 



' 1 cm 2 . Hence, for example, 
10 13 cm" 2 , 



Another explanation was proposed by lNeufeld fc Yuan! (120081 ). who encountered a simi- 
lar disagreement between mid- and near-infrared line intensities in their observations for the 
SNR IC443. They analyzed the four Spitzer IRAC band data (3.6, 4.5, 5.8, and 8.0 /xm) 
with a power-law admixture model, and found that the derived n(H 2 ) is ~ 10 6 cm -3 when 
the 3.6 /im band data are not used, while it is ~ 10 7 cm -3 when the 3.6 /zm band data are 
used. The 3.6 /im band is dominated by the v — 1 — > rovibrational lines, while other bands 
are dominated by pure rotational lines. Taken together, this implies a density that is higher 
when derived from rovibrational lines than when it is derived from pure rotational lines. 



Neufeld fc Yuan! (120081 ). however, argued that this disagreement may be caused by the 



omission of collisions with hydrogen atoms in the power-law admixture model. They noted 
that the cross section for excitation by H is several orders of magnitude greater for rovi- 
brational transitions th an it is for pure rotational transitions (see Table 1 and Figure 1 in 
Le Bourlot et al.lll999l ). Hence, with only a small fraction of H, n(H)/ri(H 2 )~ 0.025, the 
rovibrational transition can be dominated by collisions with H, rather than with H 2 , in 
the temperature range 300-4000 K. Indeed, such a fraction of n(H)/n(H 2 )~ 0.025, which 
corresponds to n(H)/n H ~ 0. 012, is expected in in terstellar clouds with n(Ho ) > 10 3 cm" 3 



(see Table 1 and Figure 1 in lSnow fc McCall 



200 61). Also. aslNeufeld fc Yuan! (120081 ) noted 



Wilgenbus et al 



some theoretical models for shock waves (e.g. 
atomic hydrogen abundances can be achieved behind shocks that are fast enough to produce 



2000j) suggest that significant 
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H2 at temperatures of a few thousand K. Since we also omitted collisions with H in the 
power-law admixture model, the higher than expected H 2 v = 1 — > S(l) intensity might 
be understood in this way. 

7. Conclusions 

We observed the shock-cloud interaction region in the SNR HB21 at near- and mid- 
infrared wavelengths, with the WIRC at the Palomar telescope and the IRC aboard the 
AKARI satellite. The IRC S7, Sll, and L15 band images and the WIRC H 2 v = 1 -> S(l) 
image reveal similar diffuse features, such as a blast wave, a bow shock, and a shocked clump. 
We chose three representative regions — Nlwake, N2clump, and N2front — and analyzed their 
infrared characteristics with several different shock models, on the basis that H2 emission 
lines are the radiation source of the shock-cloud interaction features. 

We found that the IRC colors are well explained by an admixture model of H 2 gas 
temperatures, whose infinitesimal column density varies as dN ~ T~ b dT. Three physical 
parameters — n(H 2 ), b, and iV(H2 ; T > 100 K) — were derived from this thermal admixture 
model (see Table [3]). The derived b value (~ 3) can be understood by a bow shock picture, 
whose shape is cycloidal rather than paraboloidal. However, the bow shock interpretation 
has a size-scale problem for the N2front, because the observations require that the bow shock 
be hidden within ~ 0.01 pc ~ 3 x 10 16 cm, to be unresolved in the WIRC H2 v = 1 — > S(l) 
image (T ~ 1.2"); this scale is smaller than the expected size of bow shocks. Instead of the 
bow shock picture, we propose a shocked clumpy ISM picture to simultaneously explain the 
obtained model parameters and the mid-infrared appearances. To confirm this picture, more 
robust theoretical studies and observations are required. We also compared the observed H 2 
v = 1 — ► S(l) intensity to the predicted intensity from the power-law admixture model. The 
observed H2i> = l— »0S(1) intensities are a factor of ~ 17 — 33 greater than the predicted 
ones. This excess might be caused by either an additional component of hot, dense H2 has 
(which has low total column density), or through the omission of collisions with hydrogen 
atoms in the power-law admixture model (which results in an under-prediction of the near-IR 
line intensity). 

This work is based on observations with AKARI, a JAXA project with the participation 
of ESA. The authors thank all the members of the AKARI project. Also, the authors 
thank the referee for all the comments which make this paper clearer. J.H.S thanks Ji- 
Hyun Kang for her helpful comments on the analysis. H.G.L acknowledges the partial 
support of the Korea Science and Engineering Foundation to the Astrophysical Research 
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Fig. 1. — Schematic description of H 2 level population diagrams. 
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Fig. 2. — Regions observed by AKARI. The IRC pointing-observation regions ("Cloud N" 
and "Cloud S") are overlaid as two boxes on the 1420 MHz radio continuum image of HB21, 
obtained by using the synthesis telescope at the Dominian Radio Astrophysical Observatory. 
The Palomar WIRC observations were performed toward a similar region. In this paper, we 
present the Cloud N data. The 1420 MHz radio continuum image is kindly provided by T. 
L. Landecker. 



Fig. 3. — The AKARI IRC and Palomar WIRC images of Cloud N, together with reference 
images. See Table [1] for the band definition of IRC images, {top-panel) IRC N3, N4, and 
S7 band images, (midd le-panel) IRC Sll, L15 , and L24 band images, (bottom-p anel) 1420 
MHz radio continuum iTatematsu et alJll99ph . 12 CO J = 2 ->■ 1 230.583 GHz JKoo et al. 
200lh . and H 2 v = 1 -> S(l) 2.122 fj,m images. The peak positions ("Nl" and "N2" 



where broad CO molecular lines were observed, are indicated with a '+' over all images (cf. 



Koo et al.ll200fl ). Bright point-sources were simply masked out, as the white circles indicate 



in the upper six panels. 



Fig. 4. — (top) The L15 band image of a clumpy feature at the west-end of the N2 cloud, 
designated as "N2clump" in section I3T21 (cf. Fig. ED . (bottom) Schematic description of the 
cycloidal bowshock shape. 



Fig. 5.— The IRC RGB image (left) and the WIRC H 2 v = 1 -> S(l) image (right) of 
Cloud N. The RGB image is composed of L15 (R), Sll (G), and S7 (B) band images, i.e. 
15 yum +11 yum + 7 /im. All colors scale linearly, and fully cover the dynamic range of the 
diffuse features. The three regions selected for study are overlaid on both images: Nlwake, 
N2clump, and N2front. The area defined by a black- or- white solid line is the source, while 
the area defined by a black- or- white dashed line is the background. Circular areas around 
possible point sources were excluded during the intensity measurement to avoid possible 
contamination. These areas are indicated as black- or- white circles with a red slash. Bright 
point-sources are masked out, and their positions are indicated by green circles with black 
shading. The area outlined with a yellow polygon excluding yellow circles with a red slash 
is the sampled area for the scatter plots between the IRC and WIRC images (cf. Fig. ED- 



Fig. 6. — Scatter plots between the IRC band images and the WIRC H 2 v = 1 — > S(l) 
image. Each point represents the mean intensity within an 8" x 8" region. The sampled area 
is outlined in Figure El with a yellow polygon excluding yellow circles with a red slash. The 
lines plotted over the lower-panels are the linear fitting lines. The units are MJy sr _1 and 
erg s -1 cm -2 sr _1 for the IRC and H 2 v = 1 — > S(l) images, respectively. 
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Fig. 7. — The IRC color-color diagram for Cloud N. The axes represent the ratio of the 
intensities in the corresponding IRC bands. The data points are shown by the three symbols 
in the legend (diamond, triangle, square). The expected colors for both isothermal ( §5.11) 
and power-law-thermal ( §5.2[) cases are indicated as open circles (o) and filled circles (•), 
respectively. OPR=3.0 is assumed for both cases. The different types of black lines connect 
points of equal n(H 2 ) and the LTE case. The grey solid lines connect points of equal power- 
law index (6) or equal temperature (T). The values for the power- law index and temperature 
are also indicated. 



Fig. 8. — The IRC color-color diagrams for Cloud N with the expected colors for isothermal 
cases of various OPRs (cf. §5.ip . The rest of the graph is the same with Figure [71 
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Fi g. 9. — The expected IRC colors for a partially dissociative J-shock model [open circles, 



Brand et al. 


1988: 
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1998; 
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1989 


) and the non-s 


Flower & Pineau des Forets 


1999) 



intensities in the corresponding IRC bands. The data points are shown by the indicated 
symbols in the legend. The connected open-circles have the same fractional CO abundance 
to iV(H 2 ), X co ; the four lines correspond to X co = 0, 1CT 6 , 1CT 5 , 10~ 4 , from left to right. 



The pressure increases from 10 cm K to 10 cm K, by a factor of 10 , along each line. 
Some points overlap near that for P = 10 3 cm -3 K. The annotations near the filled circles 
indicates the time passed in the non-stationary shock model. 
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Table 1. Summary of the AKARI IRC Observations 



Channel 


Filter 


Wavelength 


Imaging 


Data ID 






coverage 3. 


Resolution (r) 




(pixel size) 




(H 


(FWHM, ") 




NIR 


N3 


2.7-3.8 


4.0 


1400728 


(1.46" x 1.46") 


N4 


3.6-5.3 


4.2 


1400728 


MIR-S 


S7 


5.9-8.4 


5.1 


1400728 


(2.34" x 2.34") 


Sll 


8.5-13.1 


4.8 


1400728 


MIR-L 


L15 


12.6-19.4 


5.7 


1400729 


(2.51" x 2.29") 


L24 


20.3-26.5 


6.8 


1400729 



a Defmed as where the r esponsivity is larger than 1/e of the peak 
for the imaging mode. See lOnaka et al.l (120071 ) . 



Table 2. Results toward Cloud N 



Region 


S7 


Sll 


L15 


S7/S11 


S11/L15 


H 2 v = 1 -> S(l) a 




(MJy sr- 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 






(erg s _1 cm -2 sr _1 ) 


Nlwake 


< 0.058 


0.21±0.04 


0.36±0.08 


< 0.28 


0.59±0.16 


< 3.1xl0~ 7 


N2clump 


0.33±0.02 


0.69±0.05 


0.57±0.08 


0.47±0.05 


1.2±0.2 


(2.8±0.2)xl0~ 6 


N2front 


0.60±0.01 


0.81±0.03 


0.47±0.05 


0.75±0.03 


1.7±0.2 


(3.8±0.1)xl0- 6 



Tixtmction-corrected intensity with iV(H)=3.5 x 10 21 cm 2 . See text for details. 



Note. — The three regions are indicated in Figure [5j The upper limits are 90% upper confidence 
limits. 
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Table 3. Derived Parameters for the Power-law Admixture Model 3, and the predicted H2 

v = 1 — > S(l) intensity 



Region 


n(H 2 ) 
(cm -3 ) 


b 


iV(H 2 ;T > 100 K) 
(cm -2 ) 


predicted H 2 u = 1 -> S(l) 
(erg s -1 cm -2 sr -1 ) 


CloudN-Nlwake 
CloudN-N2clump 
CloudN-N2front 


< 5.2xl0 2 
(5.2™) xlO 2 
(1.8ig;?)xl0 3 


< 3.8 
2 q+o.3 

Z.J_q 2 


(3.6^) xlO 20 
(2.5^; 3 )xl0 20 


(8.6^) xlO -8 
(2.3lH)xlO- 7 


a See §5^1 for the 


model description 


and Fij 


$ure[7]for graphical representation, b is the power-law 



index in the equation of diV ~ T b cff. 



Table 4. Derived Contribution of H 2 line emission to the IRC bands 



Transition 


Wavelength 


Upper State Energy 


IRC 


Weight a 


% Contribution 




(fiUl) 


(K) 






N2clump 


N2front 


H 2 v = - 5(6) 


6.109 


5830 


S7 


0.346 


4 


6 


H 2 « = - S(5) 


6.909 


4586 


S7 


0.530 


45 


52 


H 2 v = - 5(4) 


8.026 


3474 


S7 


0.961 


51 


42 


H 2 w = - 5(3) 


9.665 


2504 


Sll 


0.921 


81 


84 


H 2 v = - 5(2) 


12.279 


1682 


Sll 


0.610 


19 


16 


H 2 « = - 5(2) 


12.279 


1682 


L15 


0.023 


1 


1 


H 2 v = -0 5(1) 


17.035 


1015 


L15 


1.330 


99 


99 



In units of 10 4 MJy sr V(erg s 1 cm 2 sr 1 ). See the text for the description. 



